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SELF-CONSISTENT  FIELD  MODEL 
FOR  THE  COMPLEX  CAVITY  GYROTRON 


I.  DITBiODUCTlON 

The  gyrotron  oscillator  is  currently  under  development  as  a  high  average 
power,  efficient  millimeter  wave  source.  The  gyrotron  derives  much  of  its 
high  power  capability  from  the  ability  to  operate  stably  in  a  single  mode  of 
an  over-moded  resonator.  However,  as  the  trend  toward  higher  output  power 
leads  to  increasingly  over-moded  cavities,  the  mode  density  Increases  and  mode 
competition  can  occur.  This  effect  can  seriously  degrade  the  efficiency  of 
operation  and  is  undesirable  for  other  reasons  as  well.  For  example,  the 
plasma  heating  application  for  magnetic  fusion  involves  highly  mode  and 
frequency  specific  transmission  systems  for  the  output  power.  Thus  mode 
selection  techniques  for  gyrotron  cavities  are  of  considerable  interest, 
particularly  with  respect  to  different  transverse  modes.  One  method  which  has 
been  successfully  demonstrated  in  several  devices  involves  a  sudden  change  or 
step  in  the  cavity  wall  radius.  The  "complex-cavity"  or  "step-cavity"  concept 
was  first  demonstrated  in  a  TEQ^/TEg4  configuration  by  Carmel  et  al  (1983), 
and  has  been  incorporated  in  a  60  GHz  "^EQ^/TEgg  device  (Felch  et  al  198U)  and 
in  a  140  GHz  TEg2/TEg2  device  developed  by  Varian.  The  60  GHz  Varian  device 
produces  200  kW  of  cw  power  and  it  and  the  1 UO  GHz  gyrotron  currently  under 
development  represent  the  current  state-of-the-art  of  high  average  power 
gyrotrons  in  the  U.S. 

The  essential  idea  of  the  complex  cavity  approach  is  that  a  particular 
pair  of  modes  having  the  same  azimuthal  index  m  can  be  tightly  coupled  to  form 
a  high  Q  mode  of  the  complete  cavity.  In  the  step  cavity  configuration 
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(CarMl  tt  al  1983)  mode  coupling  occurs  at  a  sudden  change  in  the  cavity  wall 
radius.  The  cavity  is  designed  such  that  the  local  waveguide  modes  of  this 
mode  pair  are  simultaneously  near  cut*off  in  the  two  cylindrical  sections  of 
the  cavity  at  the  desired  operating  frequency.  Mode  pairs  with  a  different 
mode  index  m'  will  not  be  simultaneously  resonant  in  the  cavity  sections  and 
will  form  modes  of  the  complete  cavity  which  are  less  favorable  for  gyrotron 
operation.  Such  modes  may  have  a  shorter  effective  interaction  length,  a 
lower  Q  factor,  or  lower  inter-mode  conversion  efficiency  at  the  step.  These 
effects  usually  lead  to  higher  starting  currents.  The  success  of  this 
approach  depends  on  several  factors  including  mode  density,  fabrication 
tolerances,  and  mode  conversion  at  the  cavity  step.  As  the  mode  density 
Increases,  more  than  one  mode  pair  can  become  simultaneously  resonant.  Step- 
cavity  gyrotrons  are  more  sensitive  to  machining  tolerances  than  conventional 
gyromonotrons  because  of  the  bimodal  resonance  condition.  Finally,  conversion 
into  unwanted  modes  at  the  step  can  degrade  the  mode  purity  and  Q  factor  of 
the  device.  Clearly,  theoretical  models  are  needed  which  can  Investgate  and 
quantify  these  effects. 

In  this  paper  a  theoretical  model  of  the  step-cavity  gyrotron  is 
presented.  The  model  Is  based  on  a  Self-Consistent  Field  theory  for  low  Q 
gyromonotrons  developed  previously  (Fliflet  et  al.  1982).  This  theory  treats 
steady-state  operation  in  a  single  cavity  mode.  The  axial  variation  of  the 
cavity  RF  fields  is  determined  self-consistently  with  the  electron  beam  motion 
by  integrating  the  wave  equation  simultaneously  with  the  electron  equations  of 
motion.  The  theory  accounts  for  important  aspects  of  gyrotron  operation  such 
as  modification  of  the  cavity  Q  and  resonant  frequency  by  beam  loading.  These 
effects  are  omitted  when  cold  cavity  eigenmodes  are  used  to  describe  the  RF 
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field  distribution  during  operation.  Calculations  for  step-cavity  gyrotrons 
have  to  date  been  based  on  the  use  of  cold  cavity  elgenmodes  and  It  is  of 
interest  to  explore  beam  loading  effects  for  such  devices. 

In  this  paper,  the  SCF  approach  for  gyromonotrons  is  extended  to  account 
for  a  sudden  change  or  step  in  the  cavity  diameter.  In  addition  to  the  beam 
loading  effects  noted  for  gyromonotrons,  the  method  predicts  the  mode 
conversion  efficiency  at  the  step.  Mode  conversion  at  the  step  is  calculated 
by  the  field  matching  technique  of  Masterman  and  Clarricoats  (1971).  To 
obtain  a  computationally  tractable  model  suitable  for  use  in  design  studies, 
the  following  simplifications  are  Introduced:  (i)  The  beam  Is  assumed  to 
interact  with  Just  one  TEj^^  waveguide  mode  in  each  cavity  section.  The  modes 
in  the  two  sections  will  differ  only  in  radial  index.  (2)  Non-resonant  modes 
which  may  be  excited  at  the  step  are  assumed  to  propagate  out  of  the  cavity 
without  reflection  by  the  cavity  structure.  This  is  consistent  with  the  fact 
that  such  modes  are  not  close  to  cutoff  and  will  be  weakly  affected  by  changes 
in  the  cavity  wall  radius.  With  these  assumptions,  SCF  calculations  for  the 
step-cavity  gyrotron  are  not  significantly  more  involved  than  for  conventional 
gyromonotrons.  These  assumptions  appear  to  be  consistent  with  the  step  cavity 
concept,  but  their  validity  must  ultimately  be  tested  by  comparison  with 
experiment  or  with  more  general  calculations. 

The  method  has  been  applied  to  the  analysis  of  a  TEq^/TEq^  step  cavity 
configuration  similar  to  the  NRL  experiment  (Carmel  et  al  1983).  The  results 
indicate  that  the  selectivity  of  the  cavity  for  the  m-0  mode  is  enhanced 
compared  to  a  conventional  TEqj,  gyromonotron.  The  TEq^  -  TEqj,  mode  conversion 
efficiency  at  the  cavity  step  can  be  high  (95%)  which  allows  good  output  mode 
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purity.  The  selectivity  of  the  cavity  is  optimum  when  the  m-2  competing  mode 
Is  not  cut-off  In  the  Input  section  since  in  this  case  the  Q  factor  for  this 
mode  is  degraded  by  diffraction  losses  at  the  input.  This  effect  was 
exploited  in  the  NRL  experiment  and  probably  contributed  signifigantly  to  the 
stability  of  the  device.  The  effect  of  having  the  m-2  mode  cut-off  in  the 
input  section  was  investigated. 

The  remainder  of  this  paper  consists  of  three  sections.  Section  II 
contains  an  outline  of  the  self-consistent  field  theory  and  describes  the 
present  application  of  Masterman  and  Clarricoat's  field  matching  technique  for 
waveguide  discontinuities.  Section  III  presents  the  results  of  calculations, 
and  Section  IV  contains  discussion  of  results  and  our  conclusions.  Except  as 
noted,  MKS  units  are  used  throughout. 

II.  THEORETICAL  MODEL 

A.  Self-Consistent  Field  Equations 

The  electron  beam  -  rf  field  dynamics  are  treated  according  to  a  system 
of  coupled,  non-linear,  ordinary  differential  equations  derived  previously 
(Fliflet  et  al  1982).  The  electrons  follow  helical  trajectories  due  to  a 
strong  applied  axial  magnetic  field  which  may  be  weakly  tapered  and  experience 
perturbing  rf  fields.  Space-charge  effects  are  neglected.  On  either  side  of 
the  cavity  step,  the  beam  la  assumed  to  interact  with  a  single  TE  circular 
waveguide  mode  of  the  form: 


•  Re{C  [k  J'Ck  r)9  +  imJ  (k  r)r/p]f(z)expCi(u»t-me)]} 
t  mnonoiinn  ramn 


where  u  is  the  wave  frequency,  J^,  is  a  Bessel  function  of  the  first  kind,  the 
prime  denotes  differentiation,  and  kj^^  -  Xj^^/r^  is  the  transverse  wave 
number.  The  normalization  constant  is  given  by 


Where  Xj^j^  is  the  nth  zero  of  and  r^  is  the  radius  of  the  cavity  wall  which 
is  assumed  to  be  weakly  irregular  except  at  the  step.  The  axial  dependence  of 
the  rf  fields  in  the  resonator  is  given  by  the  complex  profile  function 


f(z)  -  f(z)expC-i^(z)] 


It  is  convenient  to  use  the  normalized  momentum  variable 


U  -  TV 


where  v  is  the  electron  velocity  and  Y  is  the  relativistic  factor  given  byi 


Y-  [1  +  u^/c^i^2 


and  c  is  the  speed  of  light.  To  obtain  a  alow  time  scale  formulation  the 
transverse  momentum  is  expressed  in  terms  of  magnitude  and  phase  variables 
and  A  according  to 


u^  ♦  iUy  -  iu^  expCi(u)t/3  -  A)] 


where  s  is  the  harmonic  number.  Consider  the  case  of  an  annular  electron  beam 
propagating  In  a  circular  waveguide  structure.  Neglecting  the  Interaction 
with  the  rf  magnetic  field  components  and  considering  the  interaction  with  a 
single  harmonic  of  the  cyclotron  frequency,  the  equations  for  steady-state 
operation  are  given  by: 


)^*^3^^mn^L^Re{f  exp(isA)) 
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(7a) 


dA  Yn  C^J  .(1<_R^)3  J.(k  r,  )Re{if  expdsA)}  y  (u/s  -  Q,/Y) 
m  -  —  nm  m— a  mn  o  —  s  mn  L  —  z 


dz  u  u^ 

Z  w 


dz 


u:  dB 


2u  B  dz 
z  z 


(7b) 


(7c) 


+  ((u/c)^  '  - 

dz^ 


-2it 
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where  z  is  the  axial  spatial  coordinate,  n-e/m  is  the  charge  to  rest  mass 
ratio  for  an  electron,  and  R^  la  the  guiding  center  radius  for  a  thin  annular 


beam.  The  local  non-relati viatic  cyclotron  frequency  is  given  by 


Q 


z 


(8) 


and  the  axial  magnetic  field  is  assumed  to  be  constant  within  the  electron 
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beam.  The  Larmor  radius  rj^  -  u^/Q^.  The  first  three  equations  are  derived 
from  the  Lorentz  force  equation,  and  the  last  equation  is  derived  from 
Maxwell's  equations.  In  Eq.(7d)  1^  is  the  beam  current  and  is  the 
electron  phase  entering  the  interaction  region.  Eq3.(7a-d)  are  solved  subject 
to  appropriate  boundary  conditions  as  discussed  below. 

It  is  readily  verified  that  the  total  power  flowing  in  the  beam  and 
interacting  mode  is  conserved  by  Eqs.(7a-d). 

B.  Field  Matching  Procedure  at  Cavity  Step 

A  computer  oriented  field  matching  technique  due  to  Masterman  and 
Clarricoats  (1971)  is  used  to  treat  mode  conversion  at  the  cavity  step.  The 
present  application  involves  the  Junction  of  two  coaxial  circular 
waveguides.  The  discontinuity  plane  aperture  cross  section  is  equal  to  the 
smaller  waveguide  cross  section.  In  this  method  the  rf  fields  in  the 
discontinuity  plane  are  expressed  in  terms  of  transverse  vector  functions. 

For  a  TE  mode  the  transverse  electric  and  magnetic  fields  in  waveguide  I  are 
in  general  given  by 

m 

K  "  ^  ^®im  -  ^rm^^n 

m 

where  the  transverse  vector  functions  are  given  by 


(9a) 

(9b) 


(10a) 


9  -  Z  X  ^  tp 

V 

h  -  -  ^^'0  (10b) 

The  scalar  function  it*  satisfies  the  wave  equation 

(V^  ♦  k^)  i(»  -  0  (11  ) 

and  the  boundary  condition 


at  the  waveguide  wall.  For  a  circular  waveguide,  the  scalar  function  is  given 
by 


ton  n  1 1 


where  -  *innl'^’'w1*  vector  functions  e  and  h  satisfy  the  orthonormality 


condition 


J  ejXhj 


Where  5..  -  0  if  IkJ  or  1  if  1-j .  The  coefficient  a.-  is  the  amplitude  of  a 


mode  incident  to  the  discontinuity  plane,  a__  is  the  amplitude  of  a  mode 


reflected  from  the  plane.  Both  propagating  and  evanescent  modes  should  be 


Included.  Similarly,  in  waveguide  II  the  fields  are  given  by 


SiSiSSy 


(15a) 
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where  e  and  h  are  vector  functions  for  waveguide  II.  To  proceed,  the 
n  n 


infinite  series  representations  of  the  transverse  fields  are  truncated  after  a 


finite  number  of  terms,  say,  p  modes  in  waveguide  I  and  p'  modes  in  waveguide 


II.  The  transverse  electric  and  magnetic  fields  must  be  continuous  across  the 


aperture.  This  condition  and  the  orthogonality  of  the  vector  functions  leads 


to  a  set  of  linear  equations  for  the  mode  amplitudes: 


a.'  +  a'  -  R  (a,  +  a  ) 
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where  the  summation  over  repeated  indicies  is  understood  and  R^j^  and  Sjjjy,  are 


given  by 
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For  circular  waveguides,  Rj^jjj  and  involve  integrals  of  Bessel  functions. 


Closed  form  expressions  for  these  integrals  can  be  obtained  (Gradshteyn  and 


Ryzhlk  1965).  In  the  present  application,  waveguide  I  corresponds  to  the 
cavity  section  for  which  the  node  M  Is  nearly  cutoff  and  Interacts  with  the 
beam;  and  waveguide  II  corresponds  to  the  cavity  section  for  which  the  node  N 
Is  nearly  cutoff  and  Interacts  with  the  beam.  Except  for  the  step,  the  cavity 
consists  of  a  slightly  Irregular  waveguide  structure  which  has  little  effect 
on  modes  which  are  far  from  cutoff.  It  is  therefore  assumed  that  any  non¬ 
resonant  propagating  modes  excited  by  mode  conversion  at  the  step  will  simply 
propagate  out  of  the  system  without  reflection  at  the  cavity  Input  or  output 
sections.  Therefore,  the  incident  mode  amplitudes  are  all  zero  except  for  aj^^ 
and  a'  . 


In  the  present  model,  the  boundary  conditions  on  the  mode  M  are  applied 
at  the  cavity  Input  and  the  self-consistent  field  equations  are  Integrated  up 
to  the  cavity  step.  Only  the  beam  Interactions  with  the  mode  M  are  Included. 
This  integration  determines  the  values  of  a^i,]  and  a^|,(  since  the  interaction 
with  both  forward  and  backward  waves  has  been  Included.  To  continue  the 
Integration  of  the  field  equations  beyond  the  step,  it  is  necessary  to 
calculate  the  amplitude  coefficients  a'j^  and  The  non-zero  value 

of  a'  results  from  the  interaction  with  the  beaun  as  well  as  from  reflection 

IN 

from  the  cavity  output  structure.  These  coefficients  as  well  as  the 
amplitudes  of  other  modes  excited  at  the  step  are  obtained  as  follows: 
Eliminating  the  coefficients  a'  from  Eqs.dSa.b)  leads  to 


2S  -a'  ■  ~  6  ,)a,  6  .)a 

mn  In  mn  n  m  m,m  Im  mn  n  m  m.m  rm 


-  0,  ra»M,  and  n»N,  Eq.(20)  can  be  rewritten 


Using  the  facts  that  a 


Eq.(21)  consists  of  p  Inhomogeneous  linear  equations  for  p  unknowns  and  can  be 

solved  using  standard  matrix  methods.  The  unknowns  are  the  p-1  amplitudes 

a  ,  m«»M,  and  a'  .  The  remaining  p'  coefficients  a  are  then  given  by 
rm  In  rn 
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To  summarize,  the  self-consistent  field  approach  for  step-cavity 
gyrotrons  involves  integrating  Eqa.(7a-d)  for  the  mode  M  up  to  the  step, 
applying  the  field  matching  procedure  to  obtain  the  amplitudes  of  modes 
excited  by  mode  conversion  at  the  step,  and  then  integrating  Eqa.(7a-d)  for 
mode  N  to  the  cavity  output  region.  The  desired  solutions  are  those  for  which 
the  mode  N  satisfies  outgoing  wave  boundary  conditions  at  the  cavity  output, 
l.e.. 


df^/dz  -  -ikjff,  (23) 

where  k^  is  the  axial  wavenumber  in  the  output  waveguide.  Within  the 
assumptions  of  the  model,  the  boundary  conditions  for  the  ether  modes  are 
automatically  satisfied  by  the  matching  procedure. 

Since  ffj  is  complex,  Eq.(23)  represents  two  conditions.  Thus  the 


solution  of  the  self-consistent  field  equations  for  the  step-cavity 


configuration  -  as  for  the  gyromonotron  -  constitutes  a  two  dimensional 
eigenvalue  problem.  If  the  mode  M  has  the  form  of  a  growing  evanescent  wave 
In  the  cavity  input  section,  then  for  fixed  beam,  magnetic  field,  and 
resonator  parameters  the  eigenfunctions  correspond  to  discrete  pairs  of  values 
for  “•  ®*"®  known  a  priori  and  must  be  found  by  a 

search  procedure. 

III.  CALCULATIONS  AND  RESULTS 

Calculators  have  been  carried  out  for  a  TEq^/TEqi,  step  cavity 
configuration.  The  cavity  design  is  shown  in  Figure  1 .  The  principal 
competing  mode  pair  is  the  TE2i/TE24  combination.  The  beam  voltage  was  taken 
to  be  70  kV  and  a  beam  pitch  ratio  of  a  -  v^/v^  -  1.5  was  used.  The  cavity 
dimensions  correspond  to  operation  at  approximately  35  GHz.  A  step-cavity 
design  of  this  type  was  successfully  demonstrated  experimentally  by  Carmel  et 
al  (1982)  at  NRL.  The  cavity  is  designed  such  that  for  the  M>0  case,  the 
resonant  mode  is  close  to  cut-off  in  both  cavity  sections,  whereas  for  the  M-2 
case,  the  resonant  frequency  la  considerably  above  cut-off  in  the  smaller 
cavity  section.  The  beam  drift  tube  of  the  experimental  cavity  was  not  cut¬ 
off  to  the  TE21  mode  which  leads  to  diffraction  losses  at  the  input  and 
degradation  of  the  Q-factor  for  the  M«2  mode.  This  effect  is  readily  achieved 
in  cavities  with  the  smaller  section  resonant  for  the  TEq.|  mode  since  TEq^  and 
TE2-]  cut-off  radii  differ  significantly.  It  la  more  difficult  to  achieve  in 
cavities  where  the  smaller  section  supports  a  higher  order  TE^^^  mode  because 
for  n  >  1  there  is  much  leas  difference  in  the  cut-off  radii  for  m-0  and  m-2 


The  convergence  of  the  field  matching  procedure  at  the  step  was  checked 
by  varying  the  number  of  modes  retained  In  the  expansions  on  each  side  of  the 
step.  As  discussed  by  Masterman  and  Clarricoats  the  number  of  modes  retained 
on  each  side  should  be  approximately  equal  to  the  ratio  of  the  corresponding 
waveguide  radii.  Adequate  convergence  was  obtained  using  four  modes  on  the 
TEq^  side  and  sixteen  modes  on  the  TEq^  side. 

The  axial  profile  of  the  RF  field  corresponding  to  the  M-0  mode  is  shown 
in  Figure  2  for  typical  current  and  magnetic  field  parameters.  The  local 
efficiency  as  a  function  of  cavity  position  is  also  shown. 

Flgtire  3  shows  the  numerically  computed  oscillation  threshold  currents 
for  the  lowest  order  M>0  and  M-2  modes  as  a  function  of  magnetic  field.  The 
optimum  efficiency  operating  points  are  also  indicated.  The  m-2  calculations 
assume  the  mode  is  cut-off  at  the  cavity  input  and  thus  are  relevant  to  higher 
order  mode  configurations  such  as  TEg^/TEg^,  etc.  The  threshold  current 
curves  are  shifted  by  about  600  Gauss.  This  shift  is  much  greater  than  that 
of  the  corresponding  threshold  current  curves  for  the  TEg^  and  TE24  modes  in  a 
gyromonotron .  Unlike  the  gyromonotron  (Carmel  et  al  1982),  the  optimum 
operating  point  of  the  step  cavity  M-0  mode  is  outside  the  M-2  threshold 
current  curve.  Thus  competition  between  M-0  and  M-2  modes  should  be 
signiflgantly  reduced.  Note,  however,  that  the  minimum  starting  current  of 
the  M-2  mode  is  not  much  higher  than  that  of  the  M-0  mode  in  the  configuration 
analyzed.  Even  though  the  TE21  node  is  not  close  to  cut-off  in  the  smaller 
cavity  section,  it  is  trapped  by  coupling  to  the  TEjig  mode  and  has  significant 
interactions  with  the  beam. 


Figure  4  shows  the  efficiency  vs.  beam  current  for  several  values  of 
magnetic  field  for  the  M«0  mode.  The  maximum  calculated  efficiency  is  about 
43)^.  This  is  somewhat  less  than  can  be  achieved  in  a  well  designed 
gyromontron  (»  50t)  for  a  ■  1 .5  and  could  probably  be  Improved  by  further 
cavity  optimization. 

Mode  purity  la  an  Important  Issue  for  step-type  cavities  since 
conversion  of  operating  modes  Into  unwanted  modes  can  occur  at  the  step.  This 
effect  can  also  reduce  the  cavity  Q-factor.  The  calculated  mode  purity  at 
optimum  efficiency  la  about  952,  thus  mode  conversion  effects  due  to  the  step 
I  are  small  In  the  present  design. 


IV.  CONCLUSIONS 

The  self-consistent  field  approach  for  gyromonotrons  has  been  extended 
to  the  analysis  of  step  cavity  configurations.  The  approach  allows  the 
numerical  calculation  of  the  large  signal  efficiency,  oscillation  threshold 
currents,  loaded  Q  factor,  and  mode  conversion  effects  at  the  step.  The 
method  of  calculation  Is  relatively  efficient  -  the  computational  effort  is 
not  greatly  Increased  over  the  gyromonotron  case  -  thus  the  method  is  suitable 
for  use  in  design  optimization  studies  of  step  cavity  devices. 


The  approach  has  been  applied  to  the  Investigation  of  the  lowest  order 
M-0  and  M-2  modes  of  a  configuration  designed  for  TEqi /TEqj^  operation.  The 
magnetic  fields  values  for  optimum  effiency  operation  are  found  to  be  more 
separated  than  in  a  gyromonotron,  thus  the  competition  between  these  modes 
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should  be  reduced.  However,  the  olnlffluo  starting  currents  for  the  M-0  and  M-2 
tDodes  are  found  to  be  comparable  In  the  case  that  the  TE21  node  Is  cut-off  In 
the  electron  beam  drift  tube.  It  Is  concluded  that  a  sufficiently  large 
amplitude  TE21  node  can  contribute  significantly  to  beam  bunching  even  when 
not  close  to  the  cut-off  condition.  This  can  occur  in  step-cavity 
configurations  when  the  lower  order  competing  mode  (TE2^ )  is  cut-off  at  the 
cavity  input  and  tightly  coupled  to  the  higher  order  mode  (TE22<).  Selectivity 
against  the  M-2  mode  is  considerably  enhanced  when  the  TE21  node  is  not  cut¬ 
off  in  the  beam  drift  tube  since  in  this  case  the  overall  cavity  Q  factor  for 
the  M-2  mode  is  significantly  reduced  and,  furthermore,  the  TE21  mode  does  not 
achieve  a  large  amplitude  from  being  trapped.  It  is  argued  that  this  should 
be  a  general  feature  of  step  cavity  designs.  To  our  knowledge  this  effect  has 
not  been  pointed  out  before. 

The  calculations  show  that  conversion  into  unwanted  modes  at  the  step  is 
small  under  optimum  efficiency  operating  conditions,  of  order  95Z. 
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